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ABSTRACT 

We investigate the spherical harmonic degree (I) dependence of the "seismic" quasi-biennial oscil- 
lation (QBO) observed in low-degree solar p-mode frequencies, using Sun-as-a-star Birmingham 
Solar Oscillations Network (BiSON) data. The amplitude of the seismic QBO is modulated by the 
11 -yr solar cycle, with the amplitude of the signal being largest at solar maximum. The amplitude 
of the signal is noticeably larger for the I = 2 and 3 modes than for the / - and 1 modes. The 
seismic QBO shows some frequency dependence but this dependence is not as strong as observed 
in the 11-yr solar cycle. These results are consistent with the seismic QBO having its origins in 
shallow layers of the interior (one possibility being the bottom of the shear layer extending 5 per 
cent below the solar surface). Under this scenario the magnetic flux responsible for the seismic 
QBO is brought to the surface (where its influence on the p modes is stronger) by buoyant flux 
from the 11-yr cycle, the strong component of which is observed at predominantly low-latitudes. 
As the 1-2 and 3 modes are much more sensitive to equatorial latitudes than the / - and 1 modes 
the influence of the 11 -yr cycle on the seismic QBO is more visible in / = 2 and 3 mode frequen- 
cies. Our results imply that close to solar maximum the main influence of the seismic QBO occurs 
at low latitudes (< 45°), which is where the strong component of the 1 1-yr solar cycle resides. To 
isolate the latitudinal dependence of the seismic QBO from the 1 1-yr solar cycle we must consider 
epochs when the 1 1-yr solar cycle is weak. However, away from solar maximum, the amplitude of 
the seismic QBO is weak making the latitudinal dependence hard to constrain. 

Key words: Sun; helioseismology. Sun: oscillations, methods; data analysis 



INTRODUCTION 



It has been known since the mid 1980s that the frequencies of the 
] Sun's resonant oscillations, known as p modes, vary throughout the 
^ solar cycle with the frequencies bei ng at their highest when the so- 
. . lar activity is at its maximum (e.g. IWoodard & Noves|[l98^ . More 
^ recently a quasi-biennial (2-year) signal has b een observed in the 
frequencies of low-degree solar p modes (e.g. iFletcher et alj|2ofol 
and references therein). In this paper we investigate the same quasi- 
biennial signal in more detail. The seismic quasi-biennial signal 
■ ■ that we observe is highly correlated with the "quasi-biennial oscil- 
lation" (QBO) associ ated in oth er proxie s of the Sun's activity (e.g . 
Kanell2005l: iLi et alj|2 006: Iva novl l2007l: IVecchio & Carbonell200^ : 



1990l: iBenevolenskavall 1 998allbl ; lTobiasll2002l : IWang & Sheelevll2003l ; 
Zaa^'ashynTetalT 2010l and references therein), some of which in- 



Zagarashvili et alJl2()10l : iBadalvan & Obridkdl201l'ir and here we re 



fer to the quasi-biennial signal observed in the p-mode frequencies 
as the "seismic" QBO. We note here that the QBO observed in 
other proxies of the Sun's activity has previously been used to re- 
fer to signals with periodicities ranging from 1-3 yrs. The actual pe- 
riod varies depending on which epoch is considered and which proxy 
for solar activity is used (e.g. Kane 2005, and references therein). 
There are many different methods by which a quasi-bienni al mod- 
ulation of the 11-yr solar cycle could be generated (e.g. iHovnd 
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volve the near- surface rotational shear layer. 

The horizontal spatial structure of p modes in the Sun can be de- 
scribed by spherical harmonics with each mode being characterized by 
its spherical harmonic degree (/) and azimuthal order (m). The spatial 
structure of the oscillations means that the sectoral \m\ = I compo- 
nents of a mode are more concentrated around the equator than the 
zonal co mponents, and this is particula rly true as / increases (see, for 
example JChristens en-Dalsgaardll2002h . The sectoral components are 
predominantly sensitive to changes in activity at lower latitudes. This 
is where strong-field structures, such as sunspots, dominate. On the 
other hand the zonal (m = 0) components show a greater relative 
sensitivity to higher latitudes, where only the weak-component flux 
is present. Although the weak-component flux does increase slightly 
with the solar activity cycle it is the substa ntial increase in the strong- 
field component at low latitudes (< 45°: Ide Toma. White & HarvevI 
I2OOCII) that is largely responsible for th e observed p-mode solar cycle 
variations ( iHowe. Komm & Hill|[l999[) . Therefore, the sectoral com- 
ponents of a p mode are more sensitive to solar cycle perturbations 
than the zonal components. The sizes of the shifts of different mode 
components scale proportionately with the corresponding spherical 
harmonic components of the observed line-of-sight surface magnet ic 
field (e.g. iHowe. Komm"& Hill 19991 : IChaplin et aDl200l I2004al l3: 
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Ijimenez-Reves et alj [20041 ; IChaplin et alj|2007bl) . By examining the 
/-dependence of tiie observed changes in frequencies throughout the 
solar cycle we can gain information on the latitudinal dependence of 
the solar magnetic field. Furthermore, the latitudinal dependence of 
the 11 - yr solar cycle can be constrained usin g the frequencies of low-/ 
modes jChaplin et alj|2007 j:IChapiii]l201 ih . 

Here we investigate whether the same constraints can be put on 
the latitudinal dependence of the seismic QBO, using low-/ modes, 
by examining the frequency shifts observed in modes of different /. 
Ivecchio & Carbone (2009.) used the green coronal emission line at 
530.3 imi observed between 1939 and 2005 to examine the spatial dis- 
tribution of the QBO. They found that the signal was strong in polar 
regions and weaker around the equator. 

We make use of Sun-as-a-star (Doppler velocity) observations, 
which are sensitive to p modes with the largest horizontal scales (or 
lowest I). Resolved helioseismic observations allow higher-/ modes to 
be observed. However, in general these modes are more localized than 
are the low-/ modes and hence sensitive to local conditions. This has 
the potential to confuse the discussion of patterns on a large scale, 
which is what interests us here. Furthermore, the method used in the 
determination of frequencies for the high-/ modes is different from 
that used at low / because of the influence of mode leakage, which 
occurs because the spatial filters are not perfect due to only half of the 
Sun being observed. Although analysis based on combined data sets is 
done , great care has to be taken to minimize systematic discrepancies 
(e.g. IChaplin et alJl2004al lbl). Therefore, in this paper, we assess what 
information concerning the location of the QBO can be obtained from 
low-/ modes alone. 

The structure of the remainder of this paper is as follows: Sec- 
tion |2] describes the data used in this paper. The /-dependence of the 
raw frequency shifts is discussed in Section[3] In Section|4]we assess 
whether the observed seismic QBO is strong enough to be significant 
when the frequencies are examined for each / individually. Then, in 
Section|5] the /-dependence and frequency dependence of the seismic 
QBO is discussed. In Section[6]we use simple models to try and con- 
strain the latitudinal dependence of the seismic QBO. Finally the main 
results are summarized in Section]?] 



2 DATA 

The Birmingham Solar-Oscillations Network (BiSON) makes Sun-as- 
a-star (unresolved) Doppler velocity observations. We have used the 
BiSON data to investigate whether a seismic QBO is present in the 
frequencies of modes with < / < 3, i.e. those modes that are promi- 
nent enough in Sun-as-a-star data to allow their frequencies to be de- 
termined with enough precision and accuracy to study solar cycle vari- 
ations. 

BiSON has now been collecting data for over 30 yrs. The time 
coverage of the early data, however, is poor compared to more recent 
data. As a consequence, here we have chosen to analyse the mode fre- 
quencies observed by BiSON from 1992 October 10 to 2011 April 7. 
To study the variation in the p-mode frequencies with time, the obser- 
vations made by BiSON were divided into both 182.5-day-long inde- 
pendent subsets and 365-day-long subsets that overlapped by 91.25 d. 
T he time series used in this paper cover a shorter epoch than the work 
of lFletcher eTaP bOIOh . The start date for this analysis, i.e. 1992 Oc- 
tober 10, was chosen because after this date the duty cycles of the 
182.5-d time series were always above 60 per cent, and the average 
duty cycle of the subsets was 80 per cent. A high duty cycle was par- 
ticularly important for determining accurately the frequencies of / = 3 
(and to a lesser extent the 1 = 2) modes. This is because the heights 
of the individual m components of these modes in a frequency-power 
spectrum constructed from Sun-as-a-star data are less prominent than 



the individual in components of the / = and 1 modes. We have also 
classified the subsets as occurring during periods of either high- or 
low-surface activity. We have defined the level of surface activity to 
be "high" when the mean 10.7 cm flux observed in the 182.5 d inde- 
pendent subsets is above 100 x 10""^ Wm"- Hz"'. 

Estimates of the mode frequencies were extracted from each sub- 
set by fitting a modified Lorentzia n model to the data us ing a standard 
likelihood maximization method jFletcheretal]|2009l) . A reference 
frequency set was determined by averaging the frequencies in sub- 
sets covering the minimum activity epoch at the boundary between 
cycle 22 and cycle 23. It should be noted that the main results of this 
paper are insensitive to the exact choice of subsets used to make the 
reference frequency set. Frequency shifts were then defined as the dif- 
ferences between frequencies given in the reference set and the fre- 
quencies of t he corresponding mode s observed at different epochs in 
other subsets jBroomhall et alj|2009l) . 

For each subset in time and each / in the range < / < 3, three 
weighted-average frequency shifts were generated, where the weights 
were determined by the formal errors on the fitted frequencies: first, 
a "total" average shift was determined by averaging the individual 
shifts of the modes over fourteen overtones (covering a frequency 
range of 1.86 - 3.72 mHz); second, a "low-frequency" average shift 
was computed by averaging over seven overtones whose frequencies 
ranged from 1.86 to 2.77 mHz; and third, a "high-frequency" average 
shift was calculated using seven overtones whose frequencies ranged 
from 2.80 to 3.72 mHz. The lower limit of this frequency range (i.e., 
1.86 mHz) was determined by how low in frequency it was possible 
to accurately fit the data before the modes were no longer prominent 
above the background noise. The upper limit on the frequency range 
(i.e., 3.72 mHz) was determined by how high in frequency the data 
could be fitted before errors on the obtained frequencies became too 
large to obtain precise values for the frequency shifts due to increasing 
line widths causing modes to overlap in frequency. 



3 DEGREE DEPENDENCE OF FIfEQUENCY SHIFTS 
OVER 11-YR SOLAR CYCLE 

Fig.[T]shows the individual-/ averaged frequency shifts observed in the 
365 d subsets that overlapped by 91.25 d. The results from the 182.5 d 
subsets and the 365 d subsets are consistent and so we have chosen to 
plot the 365 d subsets only but we note that the 182.5 d independent 
subsets were used to assess the statistical significance of any signals 
in the frequency shifts as the overlaps in the 365 d subsets introduce 
artifacts in the periodicity measures we used for this purpose. The 
shaded regions in Fig.[T]differentiate between times of high- (shaded) 
and low-surface activity (unshaded). 

The 11-yr solar cycle is clearly visible in the frequency shifts. 
For each / the high-frequency range shows more sensitivity to the 11- 
yr cycle than the low-frequency range. This is consistent with previous 
observations (e.g. Broomhall et al. 2009, and references therein). 

The higher-/ modes are more sensitive to the 11-yr solar cy- 
cle than the lower-/ modes. For example, the 1 = 2 and 3 modes 
show a larger frequency shift at times of high activity than the 
1 = and 1 modes. This is in agreement with previous ob- 
servations and, as we will now explain, there are two ma in rea- 



sons for the observed / dependence (e.g. Libbrecht & WoodaK 
1990l ; IChaDlinetai]|l998l ; Ijimenez-Reves et ali,2001. ; .Salabert et al 



200% . Firstly, at fixed f requency, the mode inertia d ecreases with 



increasing de gree (e.g. i Howe. Komm & Hill 'l999|; IChaplin et al] 
,2001; Rabello-So areSjJCor zennik & Schou 2008). The mode inertia 
jChristensen-Dalsgaard & Berthomiej|l99lh is a measure of the in- 
terior mass affected by a given mode and so, as the mode inertia 
decreases a smaller volume of the interior is associated with the 
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Figure 1. Individual-/ average frequency shifts of "Sun-as-a-star" modes observed in 365 d subsets that overlapped by 91.25 d, with frequencies between 1.86 and 
3.72 mHz (total-frequency band, blue solid line, and diamond symbols); 1.86 and 2.77 mHz (low-frequency band, black dotted line, and cross symbols); and 2.80 and 
3.72 mHz (high-frequency band, red dashed line, and triangle symbols). In each panel the shaded regions denote times of high-surface activity, while the unshaded 
regions denote times of low-surface activity. 



motions generated by the mode. As such the frequencies of high-/ 
modes are more susceptible to change throughout the solar cycle than 
low-/ modes. However, the mode inertia varies very little across the 
modes examined here (at 3000 yuHz the inertia of an / = mode 
is approximately 0.3 per cent greater than the inertia of an / = 3 
mode). The second reason for the /-dependence of the solar cycle 
shifts, which is more important for low-/ modes than the mode in- 
ertia, is the latitudinal distribution of the solar magnetic field at t he 
surface (e.g.'Howe, Komm & Hiiri999';'Chaplin et alj2003ll2004al 5 
[limenez-Reves et al. 2004; Chaplin et al. 2007q). Sun-as-a-star ob- 
servations are most sensitive to the sectoral components of a mode, 
which, because of the spatial structure of the mode, show a greater 
sensitivity to equatorial regions as / increases. Therefore, since at so- 
lar maximum the solar activity is also concentrated at low latitudes, 
1 = 2 and 3 modes experience a larger shift than the / = and 1 modes 
(when observed in Sun-as-a-star data). 

Fig.[T]shows that, for each /, there is shorter-term, quasi-biennial 
structure on top of the dominant 11 -yr trend and we now study this 
structure in more detail. 



4 ASSESSING THE SIGNIFICANCE OF THE SEISMIC QBO 

We begin by assessing whether the quasi-biennial structure observed 
in the frequency shifts represent a significant QBO. The significance 
of the seismic QBO was assessed by computing periodograms of the 



raw frequency shifts. When assessing the statistical significance of the 
seismic QBO we used the frequency shifts observed in the 182.5 d in- 
dependent time series as artifacts can be introduced when using over- 
lapping subsets. Overlaps in the subsets mean that the periodograms 
are modulated by a sine squared function whose first zero was at 
1.0 yr"' i.e. the inverse of the time difference between independent 
subsets. This makes it difficult to determine whether any significant 
periodicities are present in the overlapping data, particularly above 
0.5 yr-'. 

The mean frequency shift was subtracted before each peri- 
odogram was calculated and the data were oversampled by a factor of 
10. The results are shown in Fig.|2] Also plotted in Fig. |2] are the 1 per 
cent false alarm significance levels dChaplin et aPbOOaT which were 
determined using Monte Carlo simulations. 100,000 noise time series 
were simulated, using a normal distribution random number genera- 
tor, to mimic those plotted in Fig. [T] The standard deviation of each 
point in the time series was taken as the \a error associated with the 
rotational splittings plotted in Fig. [T] Periodograms of the time series 
were calculated and the distribution of the amplitudes observed in the 
simulated periodograms was used to define the 1 per cent false alarm 
significance levels. The large peaks in each panel of Fig.[2]at 0.09 yr"' 
are the signal from the 11 -yr cycle. 

The / = 2 and 3 high-frequency-range periodograms show a 
significant peak (at a 1 per cent false alarm level) at a frequency 
of ~ 0.5 yr"'. This is t he same frequency a s the seismic QBO ob- 
served by, for example, iFletcher et al] (l2010l) . However, we note that 
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iFletcher et alj also observed a significant peak at ~ 0.5 yr"' in low- 
frequency-range frequency shifts that were averaged over < / < 2. 
However, there is no evidence for a significant signal at the same 
frequency in the low-frequency range observed here. We have im- 
proved our analysis techniques since Fletcher et al.. However, it is 
possible that this has introduced some destructive interference with 
noise that has reduced the amplitude of the signal. Alternatively the 
improved analysis could have removed constructive interference with 
noise that was artificially enhancing the observed signal observed by 
iFletcher et all Furthermore, Fletcher et al.. examined a slightly longer 
time span, that included two solar maxima, which is when the seismic 
QBO is most prominent (this will be shown in Section|5]l. 

A significant peak is observed at a marginally higher frequency 
(~ 0.55 yr"') in the high-frequency-range / = periodogram. This 
difference would be equivalent to the resolution in a non-oversampled 
periodogram and so may not be significant. The Z = 1 high-frequency 
periodogram has a significant peak at a lower frequency (~ 0.35 yr"'). 
A significant peak at the same frequency is also present in the / = 3 
low-frequency-range periodogram. We recall that the seismic QBO 
has been used to refer to oscillations with periodicities ranging from 
1-3 yrs. It is, therefore, possible that, despite the different periodicities, 
the signals observed in the different / modes can be attributed to the 
same source. 

There are two other significant peaks in the / = 3 low-frequency- 
range periodogram, at 0.60 and 0.75 yr"' respectively. The origin of 
these peaks is uncertain. However, we note that the peak at 0.75 yr"' is 
of a similar periodicity to the 1.3-yr s ignal identified in other studies 
(e.g. iHowe et alj 200( : Jimenez-Reves et al. 2003; Wang &^heeley, 
l2003( : lBroomhall et aljbol ih . 



As a check of our methodology we ran the same tests on some 
artificial data. The data were simul ated for the solar Fi tting at Low 
Angular degree Group (solarFLAG; Chaplin et 'ZI I2OO6I) . These data 
were made in the time domain and were designed to mimic Sun-as- 
a-star observations. The simulated oscillations were stochastically ex- 
cited and damped with lifetimes analogous to those of real solar oscil- 
lations. The input mode frequencies were the same in each simulated 
data set and were held stationary throughout each simulated time se- 
ries. However, the data simulated the stochastic nature of the mode 
excitation giving access to different mode and noise realizations. We 
have created a time series whose length was chosen to contain the 
same number of 182.5 d subsets as the BiSON data analyzed here. Al- 
though not shown here only one significant peak was observed in the 
periodograms of the frequency shifts. Furthermore, the peak was ob- 
served in the / = 3 frequency shifts and, as we shall explain in Section 
15.11 estimates of the Z = 3 mode frequencies are more noisy than the 
frequencies estimated for the other I. This indicates that we must be 
careful concerning the significance of the peaks observed in the I = 3 
periodograms as it is possible that the uncertainties associated with the 
/ = 3 frequencies underestimate the scatter in the estimated frequen- 



5 THE SEISMIC QBO 

We now look at the /-dependence of the QBO in more detail. In order 
to extract the seismic QBO, we subtracted a smooth trend from the av- 
erage shifts by applying a boxcar filter of width 2.5 yr. This removed 
the dominant 11 -yr signal of the solar cycle. Note that, although the 
width of this boxcar is only slightly longer than the periodicity we are 
examining, wider filters produce similar results. The resulting residu- 
als, which can be seen in Fig.[3](from the overlapping 365-d subsets), 
show a periodicity on a timescale of about 2 yrs i.e. the seismic QBO. 

For each /, the amplitude of the seismic QBO in the total- and 
high-frequency ranges varies with time. However, we note this is only 



marginally true for the 1=1 total-frequency range. Furthermore, the 
amplitude of the signal in the / = 1 total-frequency range appears to 
be smaller than the amplitude of the signal in the / = 1 low- and high- 
frequency range because, at times, the signal in the low- and high- 
frequency ranges destructively interfere with each other. It could be 
argued that the amplitude of the seismic QBO in the total-frequency 
range is also smaller than in the low- and high-frequency ranges for 
the other /, although to a lesser extent. The effect is particularly evident 
at solar minimum. 

The amplitude of the seismic QBO is larger at times of high solar 
surface activity than at times of low surface activity i.e. there is evi- 
dence for an 11 -yr envelope on top of the seismic QBO, which is in 
agreement with the results of Fletcher et al. (2010). The amplitude of 
the envelope shows some degree dependence, being particularly visi- 
ble in the / = 2 and 3 high-frequency range residuals and less visible 
in the / = and 1 residuals. Notice that, for each /, the same envelope 
is barely, if at all, visible in the low-frequency-range residuals (see the 
left-hand panels of Fig. [3]l i.e. the amplitude of the seismic QBO re- 
mains approximately constant with time. This indicates that the elfect 
of the 11-yr cycle on the seismic QBO is, to some extent, frequency 
dependent. 

The effect of the 11-yr envelope is further demonstrated by Ta- 
ble[T] which shows the maximum absolute deviations of the frequency 
residuals. For the high-frequency-range 1 = 2 and 3 modes the max- 
imum absolute deviation is greater at times of high activity than at 
times of low activity. However, the maximum absolute deviations ob- 
served in the high-frequency-range / = and 1 modes are approxi- 
mately constant. 

A similar dependence on the 11-yr solar cycle is also ev- 
ident in the QBO in other proxies of the Sun's activity (e.g. 



20 K : iMursul a. Ziege r & Vilppold 
2008: IVecch io & Carbone 2008 



iBenevo lenska^ 199^ Hathawa 

2003; Valdes- Galicia & Velasc^ . . , , ,. . , 

12009 ; Zaaarashvih e t alj|2010l and references therein). The top panel 
of Fig. |4] shows the variation throughout the last solar cycle of the 
10.7-cm radio flux (F107) emitted from the Sun and the Interna- 
tional Sun spot numb er (ISN). Man y authors (e.g. lElsworthetalJl990l ; 
IChaplin et al. 2007b: lHowell2008l . and references therein) have deter- 
mined and commented upon the good correlations between shifts in 
the p-mode frequencies and activity proxies such as the F10.7 and ISN. 
The bottom panel of Fig. |4] shows the residuals of the proxies once 
the dominant 1 1-yr signal has been removed (obtained by subtracting 
a smooth trend from the proxies). The proxies show clear evidence of 
the QBO, predominantly at times of high solar activity. 

It is, therefore, reasonable to ask: Is there any evidence for the 
seismic QBO during times of low activity? Although not shown here, 
periodograms of the high-frequency-range frequency shifts observed 
during the most recent solar minimum (the unshaded region after 2005 
in Fig. [TJ indicate that that seismic QBO is still present, with a false 
alarm significance of less than 10 per cent in the / = 0, 1, and 2 p 
modes. The seismic QBO is not significant at a 10 per cent threshold 
level in the / = 3 modes because the errors associated with the fre- 
quency shifts are larger. 

These results are not inconsistent with the conjecture that the 
magnetic flux responsible for the seismic QBO is situated in the near 
surface shear layer. Under such a scenario, when the 11-yr cycle is 
in a strong phase, buoyant magnetic flux sent upward from the base 
of the envelope could help to nudge flux produced in the near sur- 
face shear layer into shallower regions thus allowing the seismic QBO 
to be detected in the frequency signal. The 11-yr magnetic activity 
which brings the flux responsible for the seismic QBO to the surface 
is concentrated at low latitudes (5° < A < 40°). When the 1 1-yr cycle 
is at or close to solar maximum, flux originating in the near-surface 
shear layer at low latitudes will be closer to the surface than the flux at 
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Figure 2. Periodograms of the individual-/ frequency siiifts obsei-ved in different frequency bands using the 182.5 d independent subsets. In each panel the black thin 
lines are the periodograms of the low-frequency-range shifts and the red thick lines represent the high-frequency-range periodograms. The 1 per cent "false alarm" 
significance levels for the respective frequency ranges are also plotted (green dashed lines). The periodograms have been normalized so that for each periodogram the 
1 per cent significance levels are unity. Note that the underlying resolution of a periodogram of the frequency shifts that had not been oversampled is 0.054 Year"' . 



higher latitudes. The influence of the flux on the p-mode frequencies 
is strongest when the flux is close to the surface. In Sun-as-a-star data 
the frequencies determined for / = 2 a nd / = 3 m odes are dominated 
by their sectoral components {\m\ = I: IChaplirTe t al. 2004a b), which 
show a greater sensitivity to low latitudes than the zonal components 
(m = 0). 

This conjecture is consistent with the observation that the seismic 
QBO is stronger in the p modes and surface activity measures around 
solar maximum. It is also understandable that the influence of the 11- 
yr signal on the flux responsible for the seismic QBO is more visible 
in / = 2 and 3 modes than in the / = and 1 modes. 

5.1 Constraints on the 1-dependence of the seismic QBO 

Monte Carlo simulations were performed to determine how con- 
strained the observed /-dependence of the seismic QBO is for the en- 
tire epoch examined in this paper. For each degree 10,000 time se- 
ries were simulated to contain a sine wave signal with a frequency of 
0.5 yr"', and a power spectral density comparable to that observed in 
the periodograms of the solar minimum frequency shifts. Noise was 
added to the signal, using a normally distributed random number gen- 
erator, where the standard deviation of each point in the time series 
was taken as the 1 a error associated with each frequency shift plot- 
ted in Fig. [T] The simulated time series were then used to create pe- 
riodograms. The simulations imply that there is more than a 90 per 



Table 1. Maximum absolute deviations in the residuals, observed in the 365 d 
overlapping subsets, for the low- and high-frequency ranges (measured in 
//Hz). The top half of the table shows the maximum absolute deviations ob- 
served during time of high solar activity i.e. before ~ 1995 and between 1999- 
2005. The bottom half of the table gives the maximum absolute deviations 
observed during times of low solar activity i.e. between 1995-1999 and 2005- 
201 1 (also see for example Fig. [3). 



/ = /= 1 1 = 2 1 = 3 

During periods of high solar activity 

Low 0.04 ±0.03 0.04 ± 0.02 0.06 ± 0.03 0.12 ±0.03 
High 0.07 ± 0.03 0.07 ± 0.04 0.14 ±0.04 0.12 ±0.05 

During periods of low solar activity 

Low 0.03 ±0.03 0.05 ±0.02 0.06 ± 0.03 0.08 ± 0.03 
High 0.05 ± 0.04 0.07 ± 0.03 0.04 ± 0.04 0.08 ± 0.05 



cent chance that the amplitude of the seismic QBO in the / = 2 and 
3 modes is greater than in the / = modes. Similarly there is more 
than an 80 per cent chance that the amplitude of the signal in the 1 = 2 
and 3 modes is greater than in the 1=1 modes. There is a 65 per cent 
chance that the amplitude of the signal in the / = 1 modes is greater 
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Figure 3. Individual-/ residuals left after the dominant 1 1-yr signal has been removed. The plotted results were obtained from the 365 d subsets that overlapped by 
91.25 d. Right column: residuals for modes with frequencies between 1.86 and 2.77 mHz (low-frequency range). Middle column: residuals for modes with frequencies 
between 1.86 and 3.72 mHz (total-frequency range). Left column: residuals for modes with frequencies between 2.80 and 3.72 mHz (high-frequency range). First 
row: residuals for / = modes. Second row: residuals for Z = 1 modes. Third row: residuals for 1 = 2 modes. Fourth row: residuals for 1 = 3 modes. In each panel the 
shaded regions denote times of "high-surface activity", while the unshaded regions denote times of "low-surface activity". 



than in the / = modes, while the amplitude of the signal in the 1 = 2 
and 3 modes is the same, to within the associated errors. 

A similar analysis can be performed for the most recent solar 
minimum. As the / = 3 signal is not significant during the low activity 
epoch we have only considered < / < 2. The simulations showed that 
the /-dependence of the seismic QBO is poorly constrained by these 
observations and so it would be difliicult to make any inferences about 
the latitudinal distribution of the source of the seismic QBO alone. 



We also examined the /-dependence during the cycle 23 solar 
maximum i.e. 1998-2005. There is more than a 99 per cent probability 
that the seismic QBO in the / = 2 p modes is greater than in both the 
1 = and 1 p modes and there is more than a 95 per cent chance that 
the seismic QBO in the / = 3 modes is greater than in the / = and 1 
modes at solar maximum. The amplitude of the signal in the 1 = 2 and 
3 modes is the same, to within the associated errors. 

Table |2] gives the ratios between the maximum absolute devia- 
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Figure 4. Top panel: Weighted mean activity proxies observed for 365 d over- 
lapping subseries. The black solid line with the cross symbols shows the F10.7, 
which is measured in units of 10""^ W m"^ Hz"' . The red dashed line with the 
plus symbols shows the ISN. The errors on the proxies are determined using the 
standard error associated with a mean. Bottom panel: Residuals of the F10.7 and 
the ISN after the dominant 1 1-yr signal was removed. The proxy residuals were 
found using 365 d overlapping subsets. The lines and symbols are the same as 
in the top panel. The F10.7 residuals have been offset by 40 X lO"^'^ W tcT^ Hz"' 
to ease viewing. 



Table 2. Ratios of the maximum absolute deviations of the / = 1,2, and 3 
modes with the maximum absolute deviations of the 1 = modes for the high- 
frequency range residuals during periods of high solar activity. 



/= 1 



/ = 2 



/ = 3 



1.0 ±0.7 2.0 ±0.7 1.7 ±0.8 



tions of the I = 1, 2, and 3 modes and the maximum absolute de- 
viations of the Z = modes for the high-frequency range residuals 
during periods of high solar activity. Even at solar maximum the con- 
straints on the /-dependence of the seismic QBO are not stringent and 
so it may be necessary to look at the frequency shifts experienced by 
higher-/ modes before any definitive conclusions can be made. 

At solar maximum the seismic QBO appears to be modulated 
by the 11-yr solar cycle, which is known to be strongest at low lati- 
tudes, and therefore has more of an influence on the higher-/ modes. 
Therefore, the /-dependence of the seismic QBO at solar maximum 
contains information on the latitudinal dependence of both the 11-yr 
solar cycle and the seismic QBO such that the latitudinal dependence 
of the seismic QBO alone cannot be extracted. To determine the lat- 



itudinal dependence of the seismic QBO we must look at solar mini- 
mum. However, as we have shown, at this time the signal is not strong 
enough to constrain the / dependence. 

The amplitude of the QBO in the low-frequency-range 1 = 3 
residuals is noticeably larger than for the other / (see Fig. [3] and Table 
[1). This could be a genuine effect: As noted in the previous section, 
there is a well-known dependence on / in the size of the 1 1-yr solar cy- 
cle frequency shifts. However, we note that the estimates of the / = 3 
mode frequencies are more noisy than the frequencies estimated for 
the other /. Furthermore, estimates of the 1 = 3 mode frequencies are 
also influenced by the nearby, much stronger, 1=1 modes. This is 
partially reflected in the size of the error bars associated with the fitted 
frequencies (see also Figs.[T]and|3j, but, as seen in Section|4]it is still 
possible that the uncertainties are underestimated. 



5.2 The frequency dependence of the seismic QBO 



Fig.[3]shows that the seismic QBO is more evident at high frequencies 
and it could be argued that the "signal" in the low-frequency range 
is just noise. This implies that the seismic QBO signal shows some 
frequency dependence. 

We have used the 182.5 d frequencies to determine the frequency 
at which the seismic QBO stops being significant. To do this the 
weighted average frequency shifts were generated for each subset in 
time, in the manner described in Section [2] However, here we aver- 
aged the frequency shifts of the modes over four overtones only. The 
lowest frequency range for which the mean shifts were calculated was 
1.86 - 2.37 mHz, i.e. the lower limit of the frequency ranges described 
in Section[2] The next frequency range was positioned so that it over- 
lapped this range by 3 overtones i.e. 1.99 - 2.50 mHz. This process was 
repeated until the upper limit on the frequency ranges described in 
Section |2] was reached. In total the mean frequency shift was deter- 
mined for 11 frequency ranges. We have considered the 1 = 2 fre- 
quency shifts only, because over the entire epoch considered here, the 
seismic QBO is strongest in the 1 = 2 frequencies (see Fig.|2). A pe- 
riodograin of the / = 2 frequency shifts was then determined in the 
manner described in Section |4] We found that the lowest frequency 
band at which the seismic QBO was still significant at a 1 per cent 
false alarm level was 2.60-3.04 mHz. 

Residuals of the / = 2 frequency shifts were then determined 
for the frequency bands in which the seismic QBO was significant. 
The maximum absolute deviations of the residuals at times of high- 
and low-surface activity are shown in Fig.|5] For comparison purposes 
the maximum absolute deviation of the raw frequency shifts was also 
calculated, as this reflects the amplitude of the seismic 11-yr solar 
cycle. Fig.|5]shows that the frequency dependence of the seismic QBO 
is weaker than the frequency dependence of the 1 1-yr solar cycle. 

Although the 11-yr solar cycle magnetic flux is believed to be 
generated at the base of the convection zone the main influence of 
the 11-yr signal on the p-mode frequencies occurs in the upper few 
100 km of the convection zone. This is above the upper turning point 
of the lowest frequency modes examined here a nd so explains t he ob - 
served frequency dependence of the 1 1-yr cycle. lFletcheretal](l2010h 
postulated that the magnetic activity responsible for the seismic QBO 
is positioned deeper below the surface than the regions responsible for 
the 11-yr signal and as such the effect of the QBO on the p-modes 
shows a weaker frequency dependence. 

Fig. |5] also shows that the frequency dependence of the seismic 
QBO is perhaps stronger at times of high-surface activity than at times 
of low-surface activity, which reflects the frequency dependence of 
the 11-yr envelope and again supports the idea that the 11-yr signal 
influences the QBO. 
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' 0, otherwise. 
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Figure 5. Maximum absolute deviations of / = 2 residuals observed in different 
frequency bands. Each frequency band covers 4 overtones and the data are 
plotted at the frequency of the lowest overlone in each range. The red squar'es 
and the dashed line show the results for times of high-suri'ace activity and the 
blue triangles and the dotted line show the results for times of low-surface 
activity. For comparison purposes the black crosses and the solid line show 
the maximum absolute deviations observed in the raw frequency shifts, which 
represents the amplitude of the seismic 1 1-yr solar cycle. 



6 CAN THE DEGREE DEPENDENCE BE MODELED IN 
TERMS OF SPATIAL DISTRIBUTION OF MAGNETIC 
FIELD? 

We now attempt to characterise the latitudinal distribution of the seis- 
mic QBO using the / dependence of the signal de s cribed in Section 
|5] together with some simple models. IChaplinetai] ( l2007ah show that 
inferences about the spatial distribution of the 11-yr solar cycle can 
be made using simple models of how the latit udinal di s tributi on of 
solar activity affects modes of different / and m. IChapliiJ (l201ll) used 
this technique to make reasonably accurate inferences concerning the 
maximum latitude at which the strong field from the 1 1-yr solar cycle 
is observed. Here we determine whether similar models can be used 
to obtain information about the spatial distribution of the source of the 
seismic QBO. 

We have produced a simple series of models to describe the latitu- 
dinal dependence of the magnetic activity res ponsible for the seismi c 
QBO. The models are based on those used bv ' Chaplin et al. I j2007al) . 
We took the amplitude of the magnetic activity to be uniform be- 
tween some lower and upper bounds in latitude, A, and null every- 
where else. We neglect any radial dependence of the sound speed to 
simplify the problem. We assumed that the expected frequency shifts 
caused by the modeled magnetic activity were proportional to (e.g. 
iMoreno-Ins ertis & Solanki 2000) 



{I -my. 

(l + my. 



f 

Jo 



|/""(cos efhid) sinMe, 



(1) 



where /'J"(cos(?) are Legendre functions, while b(6] describes the 
magnetic activity and, more specifically, its distribution in colatitude, 
9 = 90 - A. We assumed that when the frequency residuals were zero 
the activity responsible for the seismic QBO was zero at all latitudes, 
and consequently b(d) = for all 6. When, b{d) i= for a chosen 
range of 9, the determined Svj,,, then correspond to the amplitude of 
the seismic QBO. 

For the chosen range of 9 the magnetic flux was set to some con- 
stant value (b(9) = constant t 0). The only parameters that were 
changed from model to model were the maximum and minimum lat- 
itudes at which activity was present, Ay„.^^ and A^^y, respectively. For 
each combination of /l„,ax and /i,ni„ we calculated a value for Svi,,,. To 
summarise: 



Note that, in the above scenario the total magnetic flux in the region of 
influence is proportional to the area of the Sun's surface between 
and /i,^i„ and therefore varied considerably over the range of computa- 
tions made. 

We have made no restricting assumptions over the latitudinal 
range over which the magnetic activity responsible for the QBO is 
present. So for each A^]„ in the range 0° < /imi„ < 90° we determined 
the frequency shift when /i,nax was varied between /imi„ < /ln,ax ^ 90° . 
The aim of this Section is then to find a combination of values for A^.^^ 
and /l,ni„ that best represents the observed results. 

The results in Fig. |3] show that the residuals are modulated by an 
11-yr solar cycle envelope. However, according to our earlier expla- 
nation of the seismic QBO, the amplitude of the signal is enhanced 
at times of high activity by the 1 1-yr cycle and the magnitude of the 
enhancement is dependent on /. Therefore, the /-dependence of the 
amplitude of the QBO occurs because of the spatial distribution of 
the 11-yr solar cycle rather than the flux responsible for the seismic 
QBO. However, during epochs of low surface activity the signal from 
the seismic QBO is not strong enough to allow constraints to be put 
on the /-dependence of the signal. Therefore, we use the ratios from 
times of high surface activity that are shown in Table |2] to constrain 
the models. We also note that the Monte Carlo simulations performed 
in Section \5A\ imply that at solar maximum the seismic QBO in the 
1 = 2 and 3 modes will be larger than in the 1=1 modes. 

In Sun-as-a-star data only modes where / + ra is even are visi- 
ble, because of the near-perpendicular inclination of the solar rotation 
axis with respect to the Earth. When the frequency of an / = 2 mode 
is obtained from Sun-as-a-star data, 90 per cent of the final frequency 
is determined by the \m\ = 2 components and only 10 per cent is de- 
termined by the m = component. Similarly when estimating 1 = 3 
frequencies 94 per cent of the final frequency is determined by the 
\ m\ = 3 components and only 6 per cent by the \m\ = 1 components 
jChaplinet all 12004 The 1 = 2 and 3 model frequencies shifts 
were, therefore, constructed taking these percentages into account. 

Ratios between the / = 1, 2, and 3 model frequency shifts and 
the / = model frequency shifts were computed i.e. {SviISvq) and the 
results of a few examples of the model series are shown in Fig.|6l The 
top left-hand panel of Fig.|6]shows that when = 0° the constraints 
on the ratios are satisfied when A^^^ < 45°. As A^^^ increases the /i„,ax 
at which the constraints are satisfied decreases. The models, therefore, 
imply that the QBO is strongest in equatorial regions. This is to be 
expected if the QBO is indeed modulated by the buoyant strong flux 
from the 11-yr solar cycle, which is mainly observed in equatorial 
regions. 



7 SUMMARY 

We have shown that the seismic QBO is present in the individual-/ 
frequency shifts. The main properties of the seismic QBO are: 

(i) The seismic QBO in the p-mode frequencies shows an 1 1-yr en- 
velope. The envelope is most visible in the / = 2 and 3 modes but is 
also present in the / = and 1 p-mode frequencies. A possible expla- 
nation for the difference in the amplitude of the envelope is that close 
to solar maximum, the buoyant strong flux from the 1 1-yr cycle brings 
the flux responsible for the seismic QBO to the surface, where its influ- 
ence on the p-mode frequencies is stronger. The buoyant strong flux 
from the 11-yr cycle is mainly observed at low latitudes, which the 
1 = 2 and 3 modes are more sensitive to than the 1 = and 1 modes. 
The /-dependence of the 11-yr envelope would, therefore, be influ- 
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Figure 6. Ratios of the I = 1,2, and 3 model frequency shifts and the 1 = model frequency shifts given by simple models of magnetic activity. In each panel the 
linestyles are as follows: the red dashed line represents 1=1; the black solid line represents 1 = 2; and the blue dot-dashed line represents / = 3. 



enced by the 11 -yr cycle and so would not tell us about the latitudinal 
distribution of the seismic QBO alone. 

(ii) To determine the latitudinal distribution of the flux responsible 
for the seismic QBO we must study the seismic QBO signal away 
from solar maximum. Although significant at a 10 per cent false alarm 
level, the amplitude of the seismic QBO at times of solar minimum is 
much smaller than at solar maximum. Therefore it was not possible to 
constrain the latitudinal dependence of the QBO alone. 

(iii) We have used simple models to determine the latitudinal dis- 
tribution of the flux responsible for the seismic QBO when the 11-yr 
cycle is at its maximum. The models imply that the /-dependence of 
the seismic QBO is best replicated when the magnetic flux is between 
±45°. This would be expected if, at solar maximum, the seismic QBO 
is modulated by the 11-yr cycle and the strong-field component of the 
1 1-yr cycle is mainly observed at low-latitudes. 

(iv) The seismic QBO is dependent on mode frequency but this fre- 
quency dependence is weaker than in the frequency shifts of the 11-yr 
solar cycle. The frequency dependence of the seismic QBO is stronger 
at solar maximum than at solar minimum. This result is consistent with 
the conjecture that the flux responsible for the seismic QBO is dragged 
to the surface by the buoyant magnetic flux of the 11-yr solar cycle at 
solar maximum. 

In terms of the latitudinal distribution of the flux responsible 
for the seismic QBO, the results of this paper are by no means 
conclusive. To determine whether the latitudinal dependence of the 
seismic QBO agrees with the latitudinal distribution proposed by 
IVecchio & Carbonel l l2009l) , where the signal is strong in polar regions 
and weaker around the equator, we would need to consider times of 
low surface activity. Extending this study to higher-/, where it is pos- 
sible to study a wider range of m components, may help resolve this 
matter. However, it is also possible that, as with the low-/ modes, the 
latitudinal dependence of the seismic QBO cannot be isolated from 



the latitudinal dependence of the 1 1-yr solar cycle. Observing the Sun 
throughout the next solar cycle and beyond will enable better compar- 
isons between the different / to be made. The recent solar minimum 
has been described as unusual because of its depth and longevity. It 
will be interesting, therefore, to compare the seismic QBO observed 
in the upcoming solar cycle with that observed in the last. 



ACKNOWLEDGEMENTS 

This paper utilizes data collected by the Birmingham Solar- 
Oscillations Network (BiSON). We thank the members of the BiSON 
team, both past and present, for their technical and analytical sup- 
port. We also thank P. Whitelock and P. Fourie at SAAO, the Carnegie 
Institution of Washington, the Australia Telescope National Facility 
(CSIRO), E.J. Rhodes (Mt. Wilson, Califoma) and members (past and 
present) of the lAC, Tenerife. BiSON is funded by the Science and 
Technology Facilities Council (STFC). AMB, WJC, and YE also ac- 
knowledge the financial support of STFC. We thank R. New, D.W. 
Hughes and S.M. Tobias for helpful comments and discussions. RS 
is grateful to S. Turck-Chieze and CEA/DAPNIA for providing the 
facilities required to continue her work. 



IfEFEIfENCES 

Badalyan O. G., Obridko V. N., 2011, NewA., 16, 357 
Benevolenskaya E. E., 1995, Sol. Phys., 161, 1 
— , 1998a, ApJL, 509, L49 
— , 1998b, Sol. Phys., 181, 479 

Broomhall A., Chaplin W. J., Elsworth Y, Fletcher S. T, New R., 
2009, ApJL, 700, LI 62 



© 0000 RAS, MNRAS 000, 000-000 



10 A.-M. Broomhall et al. 



Broomhall A. et al., 2011, Journal of Physics Conference Series, 271, 
012025 

Chaplin W. J., 201 1, Lecture Notes of the 22nd Canary Island Winter 
School 

Chaplin W. J. et al., 2006, MNRAS, 369, 985 

Chaplin W. J., Appourchaux T., Elsworth Y., Isaak G. R., Miller 

B. A., New R., 2004a, A&A, 424, 713 
Chaplin W. J., Appourchaux T., Elsworth Y, Isaak G. R., Miller 

B. A., New R., Toutain T, 2004b, A&A, 416, 341 
Chaplin W. J., Appourchaux T., Elsworth Y, Isaak G. R., New R., 

2001, MNRAS, 324, 910 
Chaplin W. J., Elsworth Y, Houdek G., New R., 2007a, MNRAS, 

377, 17 

Chaplin W. J., Elsworth Y, Isaak G. R., Lines R., McLeod C. R, 

Miller B. A., New R., 1998, MNRAS, 300, 1077 
ChapUn W. J., Elsworth Y, Isaak G. R., Marchenkov K. I., MiUer 

B. A., New R., Pinter B., Appourchaux T., 2002, MNRAS, 336, 

979 

Chaplin W. J., Elsworth Y, Isaak G. R., Miller B. A., New R., 2004c, 

MNRAS, 352, 1102 
Chaplin W. J., Elsworth Y, Isaak G. R., Miller B. A., New R., Thiery 

S., Boumier R, Gabriel A. H., 2003, MNRAS, 343, 343 
Chaplin W. J., Elsworth Y, Miller B. A., Verner G. A., New R., 

2007b, ApJ, 659, 1749 
Christensen-Dalsgaard J., 2002, Reviews of Modem Physics, 74, 

1073 

Christensen-Dalsgaard J., Berthomieu G., 1991, Theory of solar os- 
cillations.. Solar Interior and Atmosphere, pp. 401-478 
de Toma G., White O. R., Harvey K. L., 2000, ApJ, 529, 1101 
Elsworth Y, Howe R., Isaak G. R., McLeod C. P, New R., 1990, 
Nature, 345, 322 

Fletcher S. T., Broomhall A., Salabert D., Basu S., Chaphn W. J., 
Elsworth Y, Garcia R. A., New R., 2010, ApJL, 718, L19 

Fletcher S. T., Chaplin W. J., Elsworth Y, New R., 2009, ApJ, 694, 
144 

Hathaway D. H., 2010, Living Reviews in Solar Physics, 7, 1 

Howe R., 2008, Advances in Space Research, 41, 846 

Howe R., Christensen-Dalsgaard J., Hill F., Komm R. W., Larsen 

R. M., Schou J., Thompson M. J., Toomre J., 2000, Science, 287, 

2456 

Howe R., Komm R., Hill R, 1999, ApJ, 524, 1084 
Hoyng R, 1990, in lAU Symposium, Vol. 138, Solar Photosphere: 
Structure, Convection, and Magnetic Fields, J. O. Stenflo, ed., pp. 

359-+ 

Ivanov E. V., 2007, Advances in Space Research, 40, 959 
Jimenez-Reyes S. J., Corbard T., Palle P. L., Roca Cortes T., Tomczyk 

S., 2001, A&A, 379,622 
Jimenez-Reyes S. J., Garcia R. A., Chaplin W. J., Korzennik S. G., 

2004, ApJL, 610, L65 
Jimenez-Reyes S. J., Garcia R. A., Jimenez A., Chaplin W. J., 2003, 

ApJ, 595, 446 
Kane R. R, 2005, Sol. Phys., 227, 155 

Li K. J., Li Q. X., Su T. W., Gao P X., 2006, Sol. Phys., 239, 493 
Libbrecht K. G., Woodard M. R, 1990, Nature, 345, 779 
Moreno-Insertis R, Solanki S. K., 2000, MNRAS, 313, 411 
Mursula K., Zieger B., Vilppola J. H., 2003, Sol. Phys., 212, 201 
Rabello-Soares M. C, Korzennik S. G., Schou J., 2008, Advances in 

Space Research, 41, 861 
Salabert D., Garcia R. A., Palle P L., Jimenez-Reyes S. J., 2009, 

A&A, 504, LI 

Tobias S. M., 2002, Astronomische Nachrichten, 323, 417 
Valdes-Galicia J. R, Velasco V. M., 2008, Advances in Space Re- 
search, 41, 297 



Vecchio A., Carbone V., 2008, ApJ, 683, 536 
—,2009, A&A, 502, 981 

Wang Y, Sheeley, Jr. N. R., 2003, ApJ, 590, 1111 
Woodard M. R, Noyes R. W., 1985, Nat., 318, 449 
Zaqarashvih T. V., Carbonell M., Oliver R., Ballester J. L., 2010, 
ApJL, 724, L95 



© 0000 RAS, MNRAS 000, 000-000 



